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Abstract 

We study the half filled Hubbard model on a hypercubic lattice in infinite 
dimensions in the presence of a staggered magnetic field. An exact Ward- 
identity between vertex functions and self-energies is derived, that holds in 
any phase without broken symmetry for all values of U. Making the reasonable 
assumptions that for small enough on-site repulsion U the high-temperature 
phase is a Fermi liquid, and that in the weak coupling regime the effective An- 
derson impurity model can be studied perturbatively, we proof that Hartree- 
Fock theory and the random-phase approximation are very accurate for small 
U, and that the system develops long-range antiferromagnetic order at a finite 
temperature. 

PACS numbers: 75.10.Jm, 71.30.+h 
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I. INTRODUCTION 



The physics of the Hubbard model in d = oo dimensions is currently investigated by sev- 
eral groups 0. The essential simplification in infinite dimensions is that the self-energy 
S(k, uj) is independent of the wave vector k [|J. At the same time, the dependence on the 
frequency uj remains non-trivial, and reflects generic features of correlated electronic sys- 
tems. Recent works have mainly focused on Fermi liquid- and Mott insulating phases JIJ^J . 
Although antiferromagnetic symmetry breaking has been briefly discussed in the literature 
|5|U, and has been studied numerically in Ref. ||, it seems that the functional-integral equa- 
tions describing phases with broken symmetry have not been investigated analytically. A 
perturbative study of antiferromagnetism in infinite dimensions has recently been published 
in Ref. 0. However, within this approach one finds at weak coupling a Neel temperature T/v 
that is almost linear in the on-site repulsion U, in contradiction with Hartree-Fock theory, 
which predicts that T/v is exponentially small. In the present work we shall use a non- 
perturbative method to study antiferromagnetic symmetry breaking in the repulsive half 
filled Hubbard model in infinite dimensions. We shall mainly focus on the weak coupling 
regime, and show that in this regime Hartree-Fock theory is essentially exact. 

It is generally accepted that, at least in d > 3, the model has for any finite positive U 
antiferromagnetic long-range order (LRO) at sufficiently low temperature. However, a rigor- 
ous proof of LRO does not exist. In particular, it is not obvious that the LRO is maintained 
for arbitrarily small U > 0, because the mean-field result for the order parameter at weak 
coupling is exponentially small, and fluctuations may be important. Here we shall consider 
Hubbard models with perfect nesting, particle-hole symmetry, and finite density of states 
at the Fermi energy. We shall proof that in d = oo there exists for sufficiently small but 
finite U > a finite temperature Tn(U) > 0, such that the susceptibility for the staggered 
magnetization is divergent for all temperatures T < Tn{U). Provided that no other suscepti- 
bility diverges at temperatures larger than T/v, our result implies the existence of long-range 
antiferromagnetic order in the weak coupling regime at low but finite temperatures. 
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One way to examine antiferromagnetic symmetry breaking is to calculate the stag- 
gered magnetization M(h) in the presence of a staggered field h, take the thermodynamic 
limit, and then let h — > 0. If the spin- rotational symmetry is spontaneously broken, 
Mo = linih^oM(h) does not vanish. Alternatively, one can calculate the staggered sus- 
ceptibility x{T, U) in a parameter regime where M = 0. Antiferromagnetic symmetry 
breaking manifests itself in a divergence of x(T, U) at the Neel temperature Tn(U). In the 
present work we shall take the latter approach and show that in infinite dimensions perfect 
nesting and particle-hole symmetry imply in the weak coupling regime that Tn{U) > for 
U>0. 

The hamiltonian of the Hubbard models under consideration is given by H = TLq + U, 
with 
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where the R-sum is over N sites of a d- dimensional hypercubic lattice, h is a staggered field, 
and II = [ir, . . . , 7r] is the antiferromagnetic ordering vector (we set the lattice spacing equal 
to unity). We have defined two-component operators = [cr|, crJ, where c^ CT creates spin- 
ex fermions at site R, and a = [a x , a y , a z ] are the Pauli matrices. We allow only hoppings 
that connect different sublattices, so that e* n ' r = — 1. This implies that the band structure 
of the non-interacting model, defined by 

ek = -E*re ik " (1.3) 

r 

satisfies the perfect nesting condition 

ek+n = -£k ■ (1-4) 



Metzner and Vollhardt [jnj first pointed out that a non-trivial limit d — >• 00 is only obtained 
if the t T are properly rescaled with inverse powers of d to compensate for the increase in the 
number of neighbors in high dimensions. The nearest neighbor hopping energy should be 



scaled as t r = tjyld. For general hoppings connecting different sublattices, we require that 
t r vanishes for large d in such a way that for d — > oo the density of states at U — has a 
finite limit p(e), 

p(e) = Jim lim — ^5(e - e k ) , (1.5) 

d^oo N—>oo /V 'T^ 
k 

where the wave-vector sum is over the first Brillouin zone. Note also that the hopping 
energies depend only in the distance between the sites, and do not break the translational 



invariance of the lattice. This is sufficient to assure that 11 



p(0) > . (1.6) 

We have included the terms proportional to the density in the definition of the interaction in 
Eq. |1.2| , because then at half filling the chemical potential is exactly zero at any temperature, 
and the Hartree correction to the self-energy vanishes. The spectrum of our model has then 
particle-hole symmetry at half filling. The purpose of this paper is to show that this property, 



together with Eqs.O] and O], are sufficient to imply that for p(0)U <C 1 our model model 



has long-range antiferromagnetic order at sufficiently low but finite temperatures. 



II. ANTIFERROMAGNETISM IN INFINITE DIMENSIONS 

In this section we shall derive an exact functional-integral equation for the self-energies 
in the presence of a staggered field. Although an equivalent equation has been written down 
in Ref. ||, we use here an unconventional basis that greatly simplifies the following analysis. 

Imposing the usual periodic boundary conditions, the free part Tio of our hamiltonian 
can be transformed into block- diagonal form via Fourier transformation, 

c k . = iV- 1/2 E^ kR CR CT • (2.1) 

R 

Conventionally, one chooses the staggered field in the ^-direction, h • a = ha z . In this case 
Tto can be written as 
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where the momentum sum is over the reduced Brillouin zone of the antiferromagnet, and 
C'ko- = [ c ko-5 c k+no-]- Note that the sign of the off-diagonal elements in the quadratic form in 
Eq. |2.2j depends on the spin projection. For a derivation of the functional integral-equation 
for the exact self-energy in d = oo this introduces unpleasant technical difficulties, because 
we have to deal with two-component operators that carry in addition a spin index. A simple 
trick to avoid this difficulty is to choose the staggered field in the x-direction [|T^], h a = ha x . 
In this case Eq.[TT] can be written as 



e k -h 
— h ek+n j 

where the two-component operators are now defined by 



(2.3) 



c kf c k+n| 



(2.4) 



Note that the are composed from operators with different spin projection, and that 
the sum in Eq. |2.3| is over the full Brillouin zone. Loosely speaking, the antiferromagnetic 
symmetry breaking is now labeled by a spin flip, so that the extra spin summation in Eq. |2.2| 
can be absorbed in the second component of C^. Of course, this is only a technical point, but 
it greatly facilitates the derivation of the functional- integral equation for the exact Greens 
function. In momentum space the interaction part of our hamiltonian can be written as 

U 



U 



N 



5*(k 1 + k 2 -k 3 -k 4 ) 



ki...k 4 



c k;iT Ck 3T 



J k!,k 3 



1 



c k 2 +n| c k 4 +ni — 2 k2 ' k4 



(2.5) 



where 5k,k' is the usual Kronecker-5, and 5*(k) = J2k ^k,K, where {K} are the vectors of the 
reciprocal lattice. In Eq.[2.5| we have shifted the momentum of the last two operators by II. 
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Obviously, this leaves 5* invariant, so that the interaction term can be expressed entirely in 
terms of the components of the operator C k defined in Eq. |2.4j . 

We now introduce the imaginary-time 2x2 matrix Greens function 



<2(k, 



T — T 



< T 



C k (r)C k ( 



> 



(2.6) 



where T denotes time ordering in imaginary time, and the time evolution and thermal 
average are determined by 7Y§ + U. The corresponding non-interacting Matsubara Greens 
function is 

-l/T 







dre^ r G (0) (k, 



T 



iu n - e k h 
h iu n — e k+n 



y 
) 



(2.7) 



where u n = itT(2n + 1). The self-energy matrix is defined as usual, E(k, iu n ) = 
G ( -°^ _1 (k, iu n ) — £Z -1 (k, iu n ). The essential simplification in d = oo is that the momen- 
tum conservation can be ignored ||, and we can replace <5*(k) — > 1/iV the mathematical 
expressions for the Feynman diagrams. The self-energy is then independent of k, and must 
be of the form 



lim E(k, iuj n ) = E(za; n ) 

d— >oo 



Particle-hole symmetry implies that 



^ E(ia> n ) T(iuj n ) ^ 



(2.8) 



E(ia> n ) = -E(-io>„) 
r(zu; n ) = r(— iw n ) . 



(2.9) 
(2.10) 



The exact on-site Greens function is then given by 



(2-11) 



where we use the abbreviation G n = G(iu n ) and E n = E(io; n ). From Eq. |2.7| it is clear that 
in general the summand on the left-hand side of Eq. |2.11| depends on e k and e k+ n, and that 



therefore the summation cannot be reduced to an integration over the density of states. 
However, if we require that the non-interacting band structure satisfies the perfect nesting 
condition, Eq. |1.4| , the summand in Eq. [2.11| depends on only, so that 



G_,n 



dep(e) 



iu n - e - S n h- r n 
h — T n iu n + e — S, 



V 1 



J 



(2.12) 



The unique signature of perfect nesting is that the energy e enters in the upper- and lower 
diagonal elements with opposite sign. 

The functional-integral equation for the exact Greens function is now derived in the 
standard way |§. One defines a variational Greens function 



Q.„ — G, n 1 + E T 



(2.13) 



and a single-site impurity action 

1 



'-'imp 



cio_ n c n 



T 



rl/T 

+ U dr 
Jo 



™t( r ) - j 



(2.14) 



where C\t) = [C\(t),C\(t)] are imaginary-time two-component Grassmann fields, with 
Matsubara components C n = T f^ T dre WnT C[j\ and n a {r) = C < J.(r)C (7 (r). The functional- 
integral equation for the self-energies is then a 2 x 2 matrix equation 

l/X?{C,Ct}exp [S imp )C n Ct 



G.T. 



T JV{C,&}ex V [-S. 
— < C n Cl >s imp 



imp l 



(2.15) 



Because Si mp depends on all {S n }, the right-hand side of Eq. |2.15| is in general a non-linear 
functional of the self-energies, while the left-hand side is a non-linear function of S n . Hence, 
Eq. |2.15| is a very complicated non-linear functional- integral equation. To calculate S n , one 
should first calculate the exact Greens function of the impurity model in Eq. |2.14| for general 
choice of the and obtain an explicit expression for the right-hand side of Eq. |2.15 



After that, one should solve the resulting non-linear integral equation. Of course, such 



a calculation can only be performed numerically. However, to examine the possibility of 
symmetry breaking, it is not necessary the explicitly solve these equations. 



III. STAGGERED SUSCEPTIBILITY AND VERTEX FUNCTION 

Suppose that G n is a solution of Eq. |2.15| . In the presence of a symmetry breaking field, 
the exact Greens function is of the form 

Ig f \ 



Q-n 



\ F n G n j 



(3.1) 



where the anomalous Greens function F n is related to the staggered magnetization M(h, T) 
of the underlying Hubbard hamiltonian via 

M(h,T) = Tj2F n = Y,<CyC n > s , mp • (3.2) 

n n 

Note that antiferromagnetic symmetry breaking in the original Hubbard model translates 
into ferromagnetic symmetry breaking in the impurity model. Of course, the concept of 
antiferromagnetism is meaningless in a world consisting of two degrees of freedom. Because 
the impurity model is essentially zero-dimensional, no spontaneous symmetry breaking can 
occur in this model. Thus, the spin-susceptibility x °f the impurity model, defined via 



,1/T ,1/T 

X(T) =T dr dr' 
Jo Jo 



,1/T r l/T 

dr 

10 Jo 

X < &(T)0*C{T)Cft{?)0*C{f) >5 lmp , (3.3) 



remains finite for all values of T and U 13 1. However, \ is not identical with the staggered 



susceptibility x °f the Hubbard model, defined by 

dM(h,T) 

* (T)= i lin o^/^ • (3 ' 4) 

This is evident from the fact that the self-energies are complicated functions of the external 
field, so that the derivative in Eq.R73] does not simply produce the correlation function in 
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Eq. |3.3| . Below we shall make the relation between \ an d X precise, and show that the self- 
consistency condition, Eq. |2.15| , assures that x can diverge at a finite temperature, while \ 
remains finite. 

To derive an expression for the staggered susceptibility of the Hubbard model, let us first 
assume that the hopping energies t r in Eq. |l . 1| are only non-vanishing for nearest neighbor 
sites. In the weak coupling regime, the generalization to arbitrary hoppings, subject to the 
restrictions mentioned earlier, is trivial and will be given shortly. Setting t r = t/y/2d for r 
connecting neighboring sites, the density of states in d = oo is f| 



p{e) = po exp 



2 2 

-7rp e 



(3.5) 



with po = p(0) = [tv27r] 1 . The integration in Eq. [2.12| can then be done analytically, and 
we find that the diagonal- and off-diagonal elements of G n are given by 



G n — —irp R r , 



np R n 



h-Y r 



a, 



(3.6) 
(3.7) 



where 



^7 



- (iu n - S n ) 2 + (h- r n ) 2 



1/2 



R n = erfc 



V^Po^n exp nplQl 



(3.8) 
(3.9) 



Here erfc [a;] is the complimentary error function, and the root in Eq. |3.8| should be taken 
such that Ref2 n > 0. The leading terms of R n for small and large po|S^n| is 



(3.10) 



1 for po|O n | <^ 1 

[7rp ^n] _1 for po| ^ | > 1 
Thus, R n acts as a high-energy cutoff for frequency summations. As long as the Matsubara 
sums are dominated by the infrared regime po|^n| <S 1, the precise form of the cutoff is 
irrelevant, so that we may set R n = 1 , keeping in mind that all frequency sums are cut off 
by an energy of the order of po . 
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It is now also clear that inclusion of hoppings between the sublattices beyond the nearest 
neighbors will lead to a different value of the cutoff po, and a different functional form of 
the function R n at high frequencies. However, the low-frequency behavior of Eqs. p.6[ - p^8| will 
be unchanged, because the form of these equations carries the unique signature of perfect 
nesting, and is in this sense universal. Hence, for p U <C 1 all calculations presented below 
hold also for arbitrary hoppings, provided we substitute the value of po and the cutoff 
function R n corresponding to the particular choice of hoppings. 

To examine the stability of the paramagnetic phase, we now calculate the staggered 
susceptibility \- Substituting Eq.|3/7| into Eq.|3]2|, differentiating with respect to h and letting 
h — > 0, we obtain an exact relation between the staggered susceptibility and the self-energy 
of the Hubbard model 

X = TY,Xn (3.11) 

n 

*n = j-im = npoRn^ 1 , (3.12) 
h-^o on U n 

where we have used the fact that in a parameter regime in which the symmetry is not 
spontaneously broken lim^ Tri = 0. It is understood that fl n and R n are now defined 
by setting h = T n = in Eqs.|3.8| and |3.9|. The vertex function A n is defined via the 



Ward-identity 



A. = l-Kmf= . (3.13) 



The divergence of the susceptibility is controlled by the low-frequency behavior of the self- 
energies. Recently it has been convincingly demonstrated that in the weak coupling regime 
the phase without broken symmetry is a Fermi liquid |1]^Q]. Assuming that this is indeed 
correct, we now show that for p Q U 1 the staggered susceptibility diverges in the limit 
T -> 0. 

Particle-hole symmetry implies S(0) = 0, see Eq. |2.9| . The infrared behavior of the 
sum in Eq. |3.11| is then determined by the finite temperature effective mass- and vertex 



renormalization factors = m/m* and Z^, 
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Z- = l-^ (3.14) 

Z^l-lim^ . (3.15) 

h^o dh 



Note that Z h 1 is precisely the zero-frequency limit of the vertex defined in Eq. [3.13| . From |2. 10 



it is clear that A(iu n ) is an even function of frequency, so that the low-frequency behavior 
of Eq. [3.11| is determined by A . We emphasize that in the strong-coupling regime poll ^ 1 



the phase without broken symmetry is most likely a Mott insulator, where Fermi liquid 
theory is not valid, and the effective mass diverges |1],[2|,|| • Therefore our proof is necessarily 
restricted to the weak coupling regime. Combining Eq. |3.14] with Eq. |3.8| , we obtain for small 
frequencies 

n n ~ZZ X \u n \ . (3.16) 



Substituting Eq.gT6|into Eq|U2|, we see that in the weak coupling regime the susceptibility 



is to leading logarithmic order given by 

X = Po^[L + 0(l)} , (3.17) 
where the dimensionless factor L is in the low temperature regime PqT <C 1 given by 

£ = *T£ T ^ = ln(p T)- 1 + O(l) . (3.18) 

n \ u n\ 

Note that R n merely provides a high-energy cutoff of the order of pg 1 to the Matsubara 
sum, the precise value of which is irrelevant. Eq. |3.17| is completely general, although all 



interesting physics is hidden in the renormalization factors Z^ and Zh- Without further 
explicit knowledge of these factors, there are a priori two possibilities: The first is that 
Z^/Zh diverges at a finite T > 0. Below we shall show that in the regime < poll <C 1 this 
is indeed the case, because Zh vanishes while Z^ remains finite. The second possibility is 
that Z u /Zh remains finite. It turns out that, at least in the weak coupling regime, this is only 
true for U = 0, where Z w = Zh = 1. In this case the divergence of the susceptibility is due to 
the divergence of L as T —>■ 0. In both cases, however, \ = oo at T = 0. Note the essential 

11 



role of particle-hole symmetry to assure that Q n vanishes in the limit ui n — > 0, so that the 
infrared cutoff for the sum in Eq. |3.18| is uj = 7rT, and no other energy scale. The only way 



in which the divergence could be avoided is via a divergence of the effective mass, so that 
= 0. This possibility can be excluded, because the weak coupling phase without broken 
symmetry is by assumption a Fermi liquid. We therefore conclude that the antiferromagnetic 
susceptibility is always infinite at T = 0. This is of course due to the nesting instability 
built into our model. Hence, the zero-temperature divergence of the staggered susceptibility 
survives in the weakly interacting theory, and leads to an antiferromagnetic ground state 
with spontaneously broken symmetry. This is the first main result of this work. Of course, 
we cannot exclude the possibility that some other susceptibility exhibits an even stronger 
divergence, but this seems to be very unlikely. 

IV. WARD-IDENTITY AND PROOF OF SYMMETRY BREAKING FOR T > 

AT WEAK COUPLING 

To proof that T^(U) > 0, we now show that vanishes for U > at a finite temperature. 
To get some intuition how this might happen, let us first calculate Z^ within Hartree-Fock 
theory. In this approximation the self-energies are independent of frequency, so that S n = 
by particle-hole symmetry, and Z w = 1. The off-diagonal self-energy is given by the first- 
order exchange diagram, 

r« = -t/T5>« , (4.i) 

n 

where F^> is obtained from our general expression in Eq. |3.7| by setting S n = and T n = 
Note that Eq. |4.1| is the usual self-consistency equation for the mean-field gap A = — . 
At low temperatures Eq.jlT] reduces to 

rw = Po UL[rW - h ] , (4.2) 



where the L is defined in Eq. |3.18| . Differentiating both sides of Eq. [4.2| with respect to h and 
letting h — > yields 

12 



Z h = l-p UL . 



(4.3) 



The mean-field estimate for the Neel temperature is then obtained from Zh = 0, which yields 



T N oc p lex P ~(PoU) 



(4.4) 



From Eqs. ^.ll| , |3.12 and O it is also obvious that the Hartree-Fock approximation for the 
self-energies leads to the random-phase approximation for the susceptibility. 

To show that Hartree-Fock theory is qualitatively correct, we now derive an integral 
equation for the vertex function A n . The crucial observation is that can be determined 
by differentiation of the self-consistency equation for the self-energies. The Hartree-Fock 
self-consistency requirement, Eq. |4.1| , is the simplest possible approximation. In infinite di- 
mensions we have an exact self-consistency equation at our disposal. The obvious procedure 
is then to differentiate the off-diagonal components of both sides of the functional-integral 
equation, Eq. |2.15| , with respect to h, and taking then the limit h — > 0. After a simple 
calculation we obtain in this way an exact integral equation for the vertex function A n 

1 



27rp R n -^ 



1 + A r . 



1 



where the kernel Xn,m is the following correlation function of the impurity model, 

1 



(4.5) 



Xr. 



J^2 



< Clp x C n C\ n a x C m > Simp 



(4.6) 



Eq. |4.5| is valid for all values of U and temperatures T > T^{U). The kernel of this integral 
equation is given by the Matsubara components Xn,m of the susceptibility of the impurity 
model. The Xn,m do not have a direct physical meaning, and should be considered as 
auxiliary quantities. The susceptibility of the impurity model, defined in Eq.[3.3|, can be 
written in a form analogous to Eqs. |3.1l| and |3.12 , 



X = T J2Xn 

n 

Xn ^ ] ~Xn,m 



J^2 



(4.7) 



(4.8) 
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Of course, we do not have an analytic expression for Xn,m, but for p$U <C 1 this kernel can 



in principle be calculated perturbatively, and is free of singularities JTI . 

To determine Zh, we multiply both sides of Eq. [4.5| by T and sum over n. By construction, 
summation of the left-hand side of Eq. |4.5| yields precisely 2%, see Eqs. |3.11| and |3.12| . Recalling 
that Z^ 1 is by definition the zero-frequency limit of the vertex A n , we obtain from Eq. f4.5| 

1 



2% = X + Z h l T XnK 



np R n Q r 



(4.9) 



where 



x _ K , Ao - An 



(4.10) 



Re-arranging terms, 4.9 can be written as 



Zh = Z - x 



r.-l 



tJ2 



A, 



-Xn ~ 2Z h x 



(4.11) 



with the numerical constant Z given by 



(4.12) 



Xm Xn 

In the non-interacting limit a simple calculation of the vertex function Xn defined in 
gives 



1°) = -2Gl°) 2 = 2tt W )2 



X 



where is the value of R„ aXU = 0, and 



(4.13) 



-Z7rp -Rn ) signu;,, 



(4.14) 



is the non-interacting Greens function, see Eq. |3.6| . Combining Eq. [4.13| with Eqs. |3.11| and 



3.12, we obtain the exact identity 



2Z hX = Tj2 — 



X R y (0) 



Substituting Eq. |4.15| into Eq. |4.11| , we finally we arrive at 



(4.15) 
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^7 ^nR-n 
X n TTPO^r 



Xn \ 



(4.16) 



This equation is the central result of this work. It has the structure of a Ward-identity, 
relating a vertex function (i.e. a derivative of a Greens function) to the Greens function 
itself. This identity is a direct consequency of the fact that in infinite dimensions the exact 
local Greens function satisfies a self-consistency condition. Note that Z^ — 1 /A appears also 
on the right hand side of Eq. [4.16| . However, Ao appears exclusively through the regularized 
combination A n = 1 — (A — A n )/A , which deviates from unity only due to the non-trivial 
frequency dependence of the self-energies, and remains manifestly finite even if A diverges. 
From Eqs. |4.10| and |2.10| it is obvious that for small U and u n we have A„ = 1 + 0(U 2 uj 2 ). 
Therefore all quantities on the right-hand side of Eq. [4.16| are free of singularities, and can 
be calculated perturbatively if poU < 1. To obtain the leading logarithmic divergence of 
the sum in Eq. [4.16| , we can safely ignore vertex corrections and set X n — 1. 

From Eq. [4.16| it is clear that U = is a singular point in parameter space, because only 
in this case the enumerator in the last term vanishes, so that Z^ = Z = 1. For any finite U, 
the singular frequency dependence in the last term in Eq. [4.16| leads in the zero-temperature 
limit to an infrared divergence. 

We emphasize that Eq. [4.16| is exact and valid for all values of U and T, as long as the 
symmetry is not spontaneously broken. Note also the formal similarity with the Hartree- 
Fock result. In fact, if the right-hand side of Eq. [4.16| is expanded to first order in p U, we 
recover Eq. [4.3| . In this approximation A n is independent of n, so that X n = 1. From Eq. |4.12 
we see that this implies also Z = 1. Straight-forward perturbation theory gives to first order 
in U 



Xn = X i : ) -UG^TY,X i ° ) + 0(U> 

m 

Using then Eq. f4.14| , we see that 



(4.17) 



r 1 



Xn 



Xn\ 



n 2 p 2 U + O(U 2 



(4.18) 
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and finally 



Z h = l-p UitT J £^ ] x[l + O{p U)) . (4.19) 



This is precisely the Hartree-Fock result, which has been derived above in a much sim- 
pler way, see Eq.fO]. However, Eq. [4.16| is completely rigorous, and contains all quantum 
fluctuations neglected in Hartree-Fock theory. 

We now use the well known properties of the impurity model to show that the Hartree- 
Fock value for the Neel temperature is asymptotically exact in the limit poll — > 0. The right- 
hand side of Eq. |4.16| depends on the self-energies S n , the vertices A n , and the susceptibilities 
Xn of the impurity model, which by virtue of Eq. |2.15| can also be considered as functionals 
of the self-energies. The crucial point is now that all these quantities can in principle be 
calculated perturbatively in powers of poll, and that the perturbation theory has a finite 
radius of convergence. This follows from the facts (i) that our model is by assumption a 
Fermi liquid, so that all interaction effects can be absorbed in finite renormalization factors. 
The second cornerstone of our proof is the fact (ii) that the Anderson impurity model can 
be treated perturbatively in the weak coupling regime . The lowest order correction to the 
impurity susceptibility is given by Eq. [4.18| . For impurity models with non-interacting Greens 
function of the form = [iuj n + iAsignu; n ] _1 Zlatic and Horvatic |L4| have proven that 



the perturbation series for the spin susceptibility is absolutely convergent for any \U\ < oo, 
and that in the weak coupling regime the first few terms of the series yield an extremely 
accurate approximation of the exact result. Although the self-consistent Greens function 
G n in Eq. |2.14| will not be of the form assumed in Ref. [13], it is extremely plausible that the 



validity of perturbation theory does not depend on the precise form of the non-interacting 
Greens function 

Eqs. |4.10| and |4.12| imply that for small enough U the constant Z has a convergent ex- 



pansion with the first two terms given by 

Z = l + Zl (p U) 2 + ... , (4.20) 
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where the constant z\ is of the order of unity. Moreover, particle-hole symmetry guarantees 



that the infrared cutoff for the sum in Eq. [4.16| is f2 = nT/Z^, so that the sum is logarithmi- 
cally divergent in the limit T — > 0. Hence, the last term in Eq. [4.16| can be made arbitrarily 
large by choosing the temperature small enough, provided that U > and the effective mass 
renormalization remains finite. Consequently, for arbitrarily small but positive U there 
exists always a temperature T/v > such that Zh vanishes. 

It follows by continuity that the Hartree-Fock result becomes asymptotically exact for 
poU — > 0, and that therefore Hartree-Fock theory correctly predicts a finite staggered mag- 
netization in the weak coupling regime. The Neel temperature Tn(U) can in principle be 
calculated by expanding S n , A n , and Xn to the desired order in p U, and finding the tem- 
perature that satisfies 

Because for small poll the solution of this equation smoothly connects with the Hartree-Fock 
result, it has, at least for small enough poll, a finite solution > 0. This completes our 
proof that the Hubbard model in infinite dimensions has indeed a finite Neel temperature 
in the weak coupling regime. Combining Eqs. |3.17| and [4.16| , we obtain the following exact 
result for the staggered susceptibility in the regime T,U -C p$ , 

p Zu\n (poT)' 1 



Xn_ \ 
^nj u=0. 



(4.21) 



X 



TTpoQr, 



Xn 



(7=0. 



(4.22) 



Note that the form of this equation is very similar to the standard result of the random- 
phase approximation. In fact, to leading order in p U Eq. |4.22| reduces precisely to the 
random-phase approximation for x- 



V. CONCLUSIONS 

The existence of an antiferromagnetic instability in the weak coupling regime of the half 
filled Hubbard models with perfect nesting is not surprising. Such an instability is predicted 
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by Hartree-Fock theory in all dimensions. While in one dimension this instability is known 
to be an artifact of mean-field theory, conventional wisdom is that at least in high enough 
dimensions Hartree-Fock theory becomes very accurate. In the present work we have used 
the machinery available in infinite dimensions to show that for poll <C 1 Hartree-Fock theory 
and the random-pahse approximation are very accurate in d — oo. We have shown that the 
particle-hole symmetry of the half filled model, together with the perfect nesting property 
of the non-interacting energy dispersion, determine the behavior of the susceptibility at low 
temperatures. By differentiation of the exact self-consistency condition available in d = oo, 
we have obtained exact Ward-identities between vertex functions and self-energies. The 



central equations, Eqs. [4.5| and |4.16| , hold for all values of T and U as long as the symmetry is 
not spontaneously broken, and are the basis for our proof of the antiferromagnetic instability 
in the weak coupling regime. To show that the Neel temperature is indeed finite, we have 
assumed (i) that the high-temperature phase without broken symmetry is a Fermi liquid. 
For small enough p Q U this assumption can be justified by recent analytical and numerical 
work Another essential ingredient in our proof is (ii) that in the weak coupling 

limit the susceptibility of the impurity model can be calculated perturbatively and reduces 
to the susceptibility of the non-interacting theory in a continuous and smooth way. Both 
assumptions, (i) and (ii), are on solid grounds and can hardly be questioned. If the 
reader is willing to accept these assumptions, then our proof can be considered as rigorous. 

We have not made any prediction about the strong coupling phase. Although our fun- 
damental equations are also valid in this case, the analysis becomes more difficult, because 
self-energies and impurity susceptibilities cannot be calculated by straight-forward expan- 
sion in powers of U. In particular, for poll ^ 1 we expect that the phase without broken 
symmetry is a Mott insulator, with diverging effective mass Thus, if the ground 

state at strong coupling is an antiferromagnet (at least for nearest neighbor hopping this 
is certainly the case), then we expect a direct transition between a Mott insulator and an 



antiferromagnet as we lower the temperature. The analysis of Eqs. |4.5| , [4.1(j| , and Eq. [4.21 



which are of course also valid at strong coupling, might reveal interesting new phenomena, 



and is left for the future. 

What is the relevance of our result to finite dimensions? For d < oo the density of states 
p(e) vanishes outside a fixed interval. However, we have seen that in the weak coupling 
regime the only parameter that determines the infrared behavior of Matsubara sums is 
p(0). Therefore, we believe that in all dimensions d > 3 the Hartree-Fock theory becomes 
asymptotically exact for poll — > 0. Obviously our result does not extrapolate to d = 1. 
But in this case the Hubbard model can be solved exactly, and we know that there is no 
spontaneous magnetization for all U, even at T = ||15|| . In d = 2 Hartree-Fock theory 
is completely incorrect at T > 0, because it predicts spontaneous symmetry breaking at 
low temperatures, although the rigorous Mermin- Wagner theorem |16j tells us that this can 



happen only at T = 0. Note also that in two dimensions p(0) = oo due to Van Hove 
singularities. Thus, at weak coupling the extrapolation of the physics in d = oo to d = 2 
is not possible. For nearest neighbor-hopping and large U, the half filled square lattice 
Hubbard model is equivalent to a two dimensional quantum Heisenberg antiferromagnet, 
which seems to be ordered at T = 0. There remains the possibility that the order in the 
ground state is destroyed in the weak coupling regime [|T7],|]1J. It can be shown [jnj that even 
for arbitrarily small U the perturbation expansion is not governed by a small parameter, and 
that, in contrast to d > 3, all "perturbative" corrections to Hartree-Fock theory in d = 2 
are of the relative order of unity. Thus, d = 2 seems to be closer to d — 1, and a simple 
Hartree-Fock description of the weak coupling regime seems not to be justified. 
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